Peltier effect in Andreev interferometers 
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The superconducting proximity effect is known to modify transport properties of hybrid normal- 
superconducting structures. In addition to changing electrical and thermal transport separately, it 
alters the thermoelectric effects. Changes to one off-diagonal element L12 of the thermoelectric ma- 
trix L have previously been studied via the thermopower, but the remaining coefficient L21 which is 
responsible for the Peltier effect has received less attention. We discuss symmetry relations between 
L21 and L12 in addition to the Onsager reciprocity, and calculate Peltier coefficients for a specific 
structure. Similarly as for the thermopower, for finite phase differences of the superconducting or- 
der parameter, the proximity effect creates a Peltier effect significantly larger than the one present 
in purely normal-metal structures. This results from the fact that a nonequilibrium supercurrent 
carries energy. 

PACS numbers: 74.25.Fy, 73.23.-b, 74.45,+c 



In large metallic structures, linear-response transport 
can be described using the thermoelectric matrix L that 
relates charge and energy currents to temperature and 
potential biases. 1 The off-diagonal coefficients describe 
coupling between heat and charge currents, and indicate 
the magnitude of the thermopower and the Peltier ef- 
fect. In many cases, these coefficients are coupled by 
Onsager's reciprocal relation L a p(B) = Lp a (—B) under 
the reversal of the magnetic field BX& 

In hybrid normal-superconducting systems (see Fig.[T]) 
the Cooper pair amplitude penetrates to the normal- 
metal parts. This makes the linear-response coefficients 
L different from their normal-state values, and allows su- 
percurrents /,s,eq to flow through the normal metal even 
at equilibrium. The charge and energy (entropy) cur- 
rents I l c and P E entering different terminals can in linear 
response be written as 
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FIG. 1: Example of a 4-probe structure considered in the 
text: 5 normal- metal wires connected to each other and to 4 
terminals, of which two are superconducting (S) and two nor- 
mal (N). We take the lengths /, cross-sectional areas A, and 
conductivities a of the wires to be l/lo — (1.5,1,1.2,1,0.8) 
and Aa/Aoao = (0.8,1,0.8,1,1). Here /o, Ao and o~o are 
some characteristic values, controlling the energy scale Et — 
HD I '(I3+U + I5) 2 of the proximity effect. The system is chosen 
so as to bring out effects that depend on the magnitude of ge- 
ometrical asymmetry. In the numerics, all wires are assumed 
to be quasi-one-dimensional, / ^> \f~A. 



in terms of the biases AVj = Vj 



AT, =T 3 -T 



and the modified response coefficients L. The proximity- 
induced changes in the conductance L % l x j2ai thermal con- 
ductance L%/f (for L\{ = Li\ = and ther- 
mopower -L\y(fL\\^^dldldldl^^dL have re- 
cently been investigated both experimentally and the- 
oretically. Behavior of the remaining off-diagonal coef- 
ficient has previously been discussed in Ref. \4 us- 
ing scattering theory, but the simulations were restricted 
to small structures — making the contribution from 
electron- hole asymmetry very large. 

In this article, we note that within reasonable ap- 
proximations, in diffusive superconducting heterostruc- 
tures Eq. (pQ) can be generalized to the non-linear regime 
by defining an energy-dependent thermoelectric matrix 
L l £p{E). We show that this quantity satisfies an Onsager 



reciprocal relation L 1 ^ 



, afj (E, B) = L 3 * a {E, -B) under the 
reversal of the magnetic field B and the phase argA of 
the superconducting order parameter, whenever i and j 



refer to normal terminals. We also show how the proxim- 
ity effect modifies L21 giving rise to a large Peltier effect 
and discuss how it could be experimentally detected. 



Qualitatively, one can understand the origin of 
proximity-induced thermoelectric effects by noting that 
charge current consists of a quasiparticle component 
and a supercurrent component. That the latter is 
strongly temperature dependent in proximity structures 
then leads to a finite L12 coefficient, 9,11 via a mech- 
anism analogous to charge imbalance generation in 
superconductors ] 18 i 19 Assuming Onsager symmetry, one 
would also expect that L21 is finite. The actual form of 
the coupling can be seen by inspecting the quasiclassi- 
cal transport equations (Eqs. (|4j) below), or by studying 
their near-equilibrium approximation in a diffusive nor- 
mal metal under the influence of a weak proximity effect 



2 



(see for example Ref. H2I ): 



equations 



V • J c = , V • J E = , 
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Here, SV and 5T are deviations of the (effective) local po- 
tential and temperature from equilibrium, and T is the 
ambient temperature. The first terms in charge and en- 
ergy current densities J c , Je can be considered the quasi- 
particle current and the rest the (non-equilibrium) su- 
percurrent; a and a t h are the proximity- modified charge 
and thermal conductivities, Js,eq is the equilibrium su- 
percurrent density, and T is a small factor associated 
with non-equilibrium super current. Although Eqs. ([2]) 
are not of the usual form of normal-state transport 
equations one can see that a variation 5T generates a 
change in the charge current, and that non-equilibrium 
(SV ^ 0) supercurrent carries energy current. The cor- 
responding response coefficients in Eqs. (|2b|) and ([2c|) 
are not independent, which is a signature of the On- 
sager symmetry. Comparing the magnitude of the co- 
efficients, it turns out that at low temperatures a large 
part of the thermoelectric coupling indeed arises from 
the temperature-dependence of Js,eq- At high tempera- 
tures where it vanishes exponentially, other sources be- 
come more important ] 1Q i 12 i 13 

However, validity of Eqs. ([2]) is somewhat restricted, 
since they are correct only in the linear response and to 
the first order in the proximity corrections, additionally 
assuming that the energy gap |A| of the nearby super- 
conductors satisfies feT <C |A|. For quantitative calcu- 
lations of the multiterminal transport coefficients, and to 
evaluate the proximity-corrected coefficients in Eq. ([2]), 
we start from the full non-equilibrium formalism. 

The superconducting proximity effect can be described 
using the quasiclassical BCS-Gor'kov theory j 2Q i 21 Here, 
we concentrate on diffusive normal-metal structures that 
are connected to superconducting and normal terminals, 
and neglect any inelastic scattering. The model then 
reduces to the Usadel equations ^ 21 ! 22 whose first part, 
the spectral equations, can in this case be written as 



DV 2 = -2i(E + i0 + ) sinh + ^ sinh(2<9) , (3a) 
V • (-v s sinh 2 0) = , v s = D(Vx ~ 2eA/h) . (3b) 



They describe the penetration of the superconducting 
pair amplitude F = e %x sinh 6 into the normal metal. 
We denote the diffusion constant of the metal here by 
D, and the magnetic vector potential by A. At clean 
contacts to bulk superconductors, the pairing angle is 
6 = artanh(|A| / E) and the phase x — arg A, where A is 
the superconducting order parameter. Transport prop- 
erties are in turn determined by kinetic Boltzmann-like 



DV ■ f T f = Kfr + £>(V ■ j s )h, DV ■ f L f = 0, (4a) 
f T f = V T Vf T + TVf L + jsh , (4b) 
f L f ee V L V/ L - 7V/ T + jsh , (4c) 

that describe the behavior of the antisymmetric and sym- 
metric parts f L (E) = f(jJL S -E)-f(ns+E) and f T (E) = 
1 — /(/is — E) — f(ns + E) of the electron distribution 
function. They are defined with respect to the potential 
of the superconductors, chosen below as us = 0. The 
spectral supercurrent js, the diffusion coefficients 
T>t, % and the condensate sink term 1Z are functionals 
of and x, having the symmetries T^l/t{x) = ^l/t(~ x)> 

HX) = -T[-X), js(x) = -Js(-X), n{ X ) = K(-x)W 
In normal metals, V • js = TZ = 0. Observable cur- 
rent densities are finally related to the spectral currents 
t L/T f through 



2e 2 L 



(5) 



and the heat current density is Jq = Je — VJ C at the 
terminals. Below, we also assume that all contacts to 
terminals are clean and of negligible resistance: in this 
case all quantities are continuous at the interfaces, except 
at superconductors for E < |A| the boundary condition 
for the kinetic L-mode is changed to n • T^f = 0, where 
n is the normal to the interface. 

It is important to note that the last two terms in 
Eqs. (|4]) mix the L and T modes and cause thermoelec- 
tric effects: near equilibrium, they lead to the coupling 
terms in Eqs. ([2]). Away from linear response, a non- 
equilibrium modification of the distribution function / 
due to the mixing 2 ^ has also been experimentally ob- 
served in Ref. Izj . 

The aim in the following is to calculate the thermo- 
electric coefficients L l K starting from Eqs. (|4]). However, 
as with the charge conductance, it is useful to first de- 
fine corresponding energy-dependent thermoelectric co- 
efficients L l ^p(yE). Since the kinetic equations are linear, 
it is possible to write the currents entering different ter- 
minals as 



/oo 
dE J2L%(E)f^E) : 
-CO a- 



(6a) 
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/CO 
d££^L^)4(E), (6b) 

where /3 G {T, L}, j runs over all terminals, and is the 
a-mode distribution function in terminal j. This spec- 
tral thermoelectric matrix L l ^{E) is the quasiclassical 
counterpart to the P-matrix in Ref. 4. More explicitly, 
L l £p{E) can be defined as the a-mode current seen in 
terminal i that a unit excitation of mode (3 in terminal j 



generates at energy E: 



Ltl 



dS n • f 0,-0- 



3,P 



(7) 



Here, S{ is the surface of terminal i and n the cor- 
responding normal vector. The two-component func- 
tion ip J P is assumed to satisfy the kinetic equations (|4]) 
with the electron distribution functions f l a in termi- 
nals replaced by S a pSij. The linear-response coefficients 
L are directly related to L(E) via Eq. (|6|), for exam- 
ple Lu = / ^EL TT (E) sech 2 (^ T ) and L 2 i = 

The spectral thermoelectric matrix depends only on 
and x? but not on the distribution functions at the ter- 
minals. Knowing the energy dependence of this matrix, 
one can directly evaluate currents also away from lin- 
ear response, if changes in the order parameter A and 
any inelastic scattering can be neglected. The matrix 
L % ^p(E) one can evaluate numerically once and x have 
been solved, and it offers a feasible way to find the re- 
sponse of the circuit to different types of excitations in 
the terminals. 

An Onsager reciprocal relation for L l ^{E) follows from 

the fact that the differential operator C in Eqs. (@J), £f = 
0, has the property 



£(B)t = (-V) 



n 
o 



^ TP T V(-v) + (-v) 
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= C(-B), 

due to the symmetries of the coefficients under reversal 
of the phases x? arg A and the magnetic field B. Below, 
whenever we discuss reversal of £?, also reversal of the 
phases is implied. Integration by parts now shows that 
for any volume Q and two-component functions (/>, i/j, we 
can write 



/ dV [^t<j>- 
Jn 



tf&il)] = [ ( 
Jon 



dS n • J , 



(9) 



where dft is the boundary of Q. For the differential oper- 
ator here, the flux J = ^Y{B)<j>- ^f(-B)ip - j s ^ ] "<7i0, 
a\ being the first spin matrix. Now, we choose Vt to be 
the whole conductor, with <j> and such that <p_ = 
satisfies the conditions in the calculation for L l ^{E^B) 

and ijj = the conditions for L J p a (E, —B). When both 
i and j refer to normal terminals, we then find 



= / dSh-J 

Jan 

= [ dSn-f a (B)0- / dSn.f>(-B)^, 
J Js, 



(10) 
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FIG. 2: (Color online) Elements for i,j = 1,2 and a, /? = 
L, T for the asymmetric interferometer in Fig. [T] in units of 
I/R2. The order parameters in the superconducting termi- 
nals have arg A 3 - arg A 4 = 0.547T and |A 3 | = |A 4 | = 40E T . 
Thick blue solid line is L 11 , dash-dotted red line L 22 , black 
dotted line L 12 and magenta dashed line L 21 . Note the sym- 
metry L^(E) = (-l) 1_<5a/3 ^(^). Approximations found 
by solving Eqs. J4| to first order in js and T are shown as 
thin solid lines — in general they are indistinguishable from 
the exact numerical results. 



using the boundary conditions imposed on <p and ip, and 
the fact that js = at normal terminals. Comparison of 
Eqs. ([TO]) and (|7J) reveals a reciprocal relation 



L^(E,B)=L^ a (E r 



-B). 



(11) 



This implies that phase differences in the order param- 
eter will be similar sources for quasiclassical Peltier and 
Thompson effects as they are for the thermopower dis- 
cussed in Refs. I^ ITollTlllT2lfT3l . Similar relations exist also 
in the scattering theory^ 

The form of Eqs. ([6]) also implies that L(E, B) has the 
symmetries 

^^L l r^ L (E) = for normal terminal z, (12a) 

3 

Y J L i l L {E)=0, (12b) 

3 

L%(E, -B) = (-iy-^L%(E, B) , (12c) 

since the charge current to any normal terminal and the 
entropy current to any terminal must vanish at equilib- 
rium for all temperatures. Equation (|12cp follows essen- 
tially from the electron-hole symmetry assumed in the 
quasiclassical theory, leading to TlJ ► T^/, Trf 
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FIG. 3: (Color online) Peltier coefficients Unn, Uns for the 
same parameters as in Fig. [2] Approximations (|14a|) . (|14b[) 
are shown as thin lines - deviation from the exact result is 
due to neglecting T. 



— f rpf under the transformations B —B, » — /t- 12 
This makes the diagonal coefficients symmetric in B and 
the off-diagonal ones antisymmetric. However, there are 
some experimental results 8,14 where the latter symmetry 
does not hold. Such observations cannot be explained 
with the quasiclassical theory applied here. 

Consider now the application of the formulation above 
in the structure in Fig. [TJ We solve the spectral equa- 
tions (|3j) in this structure numerically and calculate the 
spectral thermoelectric matrix from the solutions. Be- 
havior of the two coefficients important for thermoelec- 
tric effects, spectral supercurrent js and the coefficient 
T, is discussed for structures of this type for example in 
Refs. E2II25I . Resulting elements of L^JE) are plotted as 
a function of E in Fig. [2] — the energy scale is given by 
the Thouless energv E T = HD/(l 3 + U + h) 2 - The diag- 
onal elements Ltt(E) and Lll(E) are spectral charge 
and energy conductances^ At E > |A|, energy current 
can enter also the superconductor, which is visible as a 
rapid change in the L /^-coefficient. The off-diagonal co- 
efficients qualitatively follow the energy dependence of 
the spectral supercurrent js which gives the most visi- 
ble contribution. Moreover, the elements of the matrix 
clearly exhibit the symmetries ([TT]) and ([T2]) . 

The finite coefficient L l [ T {E) leads to a Peltier effect: 
assume that the terminals are at a constant temperature 
T l = T and biased at potentials chosen so that a current 
I c flows between terminals 1 and 2, l\ — — I 2 = I c . Then, 
the Peltier linear-response coefficient for this system is 
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dll 



dl\ = LWWr - L\ 2 W 2 
dll LWWx - L\\W 2 



(13) 



where Wj = (L^+L^)- 1 . 



We can also define the Peltier 
coefficient Uns = ^dl\jdl\ corresponding to the current 
configuration l\ — I 2 = I c /2. 

The magnitude and temperature dependence of II is 
shown in Fig. [3j For a typical Thouless energy of Et = 
200mK&£ of an Andreev interferometer, the Peltier co- 
efficients would be I II 7v a/" I ~ lOOnV and III/vsl ~ 1 /iV 
at T ~ 200 mK. For comparison, Peltier coefficients for 
purely normal- metal junctions at these temperatures are 
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FIG. 4: (Color online) Temperature of terminal 1 in Fig. [T] 
(a) for ip = argA3 — argA4 = ±0.547r (dashed and solid 
lines) and current configurations corresponding to II n n and 
Uns (red and blue lines), as a function of I c . (b) T\ for 
eRiIc/Er — 0.15, as a function of ip. The results are 
calculated assuming other terminals are at the temperature 
T — 2Er/kB- Deviation of T\ from T originates from Joule 
heating and the oscillation of the proximity- Peltier effect. 



of the order II = T(S B - S A ) - 0.2 K x 10 *v/ K = 2 nV. 
The interferometer induces a significantly larger II. 

The above Peltier effect is related to the thermopower 
discussed in Refs. l9lllll. We indeed find the Kelvin rela- 
tions Unn — TS NNi TIns — TSns, which follow from 
the Onsager symmetry. Similarly as in Ref. [l2|, within 
the assumptions where Eqs. ([2j) apply, one can also find 
simple approximations up to first order in js'- 
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(14b) 



Here, Is,eq = 4f dE js tanh 2fc ^ T is the equilibrium 
supercurrent. The above also shows the dependence on 
the asymmetry for Unn and the proportionality to the 
supercurrent — for this contribution to the effect. 

Finite Peltier coefficients allow for cooling one of the 
terminals by driving electric current. Assume the ter- 
minal is small enough, such that the power flowing 
into the phonons is small compared to the heat cur- 
rent carried by electrons. The temperature change 
is then limited by the Joule heat generated in the 
wires: the heat current is Iq = — G^ATi — 2UnsIc + 
e(/^) 2 /G, G and G t h being electrical and heat con- 
ductances. The maximum cooling effect then is, in a 
rough estimate assuming that Wiedemann- Franz law ap- 
plies, ATi = -(3/7r 2 )(e 2 U 2 NS /k B T)/k B 0.3 mK for 

Et = 200mKk B . Numerical calculation in the struc- 
ture of Fig. [1] yields cooling AT ~ — 0.4mK, as shown in 
Fig.H 

One point to note is that also the ^-symmetric oscil- 



5 



lation of the thermal conductance 5,6 contributes to the 
temperature change, although this is significant only at 
temperatures small compared to EtI^b- In the absence 
of the Peltier effect, AT would hence be symmetric in 
B and always positive. The proximity-Peltier effect al- 
lows negative temperature changes and also breaks the 
symmetry, which makes the antisymmetric part T\(B) — 
T\(—B) the experimentally interesting signal. In the 
structure of Fig. [T] the oscillation amplitude can be of 
the order of lmK for E T = 200mK/c£>. (See Fig. g)b.) 
Temperature changes of this order can be experimentally 
resolved in mesoscopic structures^ so that the detection 
of the effect simply via observing AT should be experi- 
mentally viable. In addition to the off-diagonal thermo- 



electric coefficients L12, T21, it would also be interesting 
to study the Onsager reciprocity for Ljl T (E) via differ- 
ential conductances in multi-terminal structures. 

In summary, we have studied charge and en- 
ergy transport and its symmetry relations in normal- 
superconducting hybrid structures. We show that a large 
Peltier effect controlled by the phase difference over a 
Josephson junction can arise, partly due to co-flowing 
quasiparticle and supercurrents. This complements pre- 
vious studies of a related effect in the thermopower. 
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